BEHAVIOR DOMINATED BY SLOW PARTICLES IN A DISORDERED 
ASYMMETRIC EXCLUSION PROCESS 

ILIE GRIGORESCU, MIN KANG, AND TIMO SEPPALAINEN 

Abstract. We study the large space and time scale behavior of a totally asymmetric, 
nearest-neighbor exclusion process in one dimension with random jump rates attached to 
the particles. When slow particles are sufficiently rare the system has a phase transition. 
At low densities there are no equilibrium distributions, and on the hydrodynamic scale 
the initial profile is transported rigidly. We elaborate this situation further by finding the 
correct order of the correction from the hydrodynamic limit, together with distributional 
bounds averaged over the disorder. We consider two settings, a macroscopically constant 
low density profile and the outflow from a large jam. 



1. Introduction 

We study a totally asymmetric, nearest-neighbor exclusion process on the one-dimensional 
integer lattice Z with random rates attached to the particles. The process is studied through 
the labeled particle configuration. The particles are indexed by integers in an increasing 
fashion. The position of particle i at time t is denoted by an integer-valued random variable 
(7i(t). The exclusion rule stipulates that ai{t) < cn+i(t) for all i S Z and all t > 0. 

At the outset each particle Oi receives its jump rate pi which then remains fixed through- 
out the dynamics. The rates p = {p{\ are i.i.d. random variables with common distribution 
F. F is supported on (c, 1] for some c > 0, and we take c to be the left endpoint of the 
support of F. In other words, F(p) = for p < c, F(p) > for p > c, and F(l) = 1. We 
also assume F(c) = so no particle has c as its intrinsic jump rate. 
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Once the rates have been fixed and an initial configuration a = {<j{ : % € Z) specified, 
the process a(t) = (cj(f) : i G Z) evolves in the usual way: each particle Oi carries its own 
Poisson clock of rate pi, and whenever the clock rings, Uj advances one step to the right 
provided the next site to the right is vacant. 

It is also useful to consider the gaps f]i(t) = Oi+\{t) — <7j(i) — 1. The process r](t) = {f]i{t) : 
i G Z) is a zero-range process with random rates attached to the spatial positions. The 
jump rule is that whenever a particle is present at position i [rji > 1], one particle is moved 
from i to i — 1 at rate pi. We can also view this system as a series of tandem queues where 
queue i is served at rate pi, and customers departing queue i immediately join queue i — 
The gap variable rji(t) is the queue length and the particle increment cxj(i) — crj(O) is the 
departure process from queue i. 

Fix the rates p = {pi}- Given any a £ [0, c], the product distribution P p with geometric 
marginals 

(1.1) pp [rji = k] = f 1 -l\fl\\ k = 0,1,2,... 

is an invariant distribution for the gap process r)(t). In this equilibrium, each particle 
motion is marginally a Poisson process with rate a. More precisely, for each i the increment 
(Ti(t) — (Ti(s) is Poisson with mean a{t — s). This is a consequence of Burke's theorem 
from queueing theory, according to which the departure process of an M/M/l queue in 
equilibrium is a Poisson process. 

When F is suitably chosen this model manifests a phase transition. Here is a way to 
approach it. Given a 6 [0, c], the (annealed) mean gap in equilibrium is 

u = [ EV[ Vi ] F® z (dp) = [ — dF{p). 
J J(c,i] P-a 

The common velocity a of the particles cannot exceed c because there are particles whose 
intrinsic rates come arbitrarily close to c. Thus the maximal mean gap u* is defined by 
letting a /* c, in other words 

(1.2) u*= [ ^—dF{p). 

J(c,i] P ~ c 

If this integral is finite, there is a critical gap size u* < oo such that the geometric product 
equilibrium distributions do not exist for mean gaps u > u* . Equivalently, there is a positive 
critical density p* = (1 + n*) _1 for the exclusion particles such that the product equilibria 
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for the gaps do not exist at low densities p < p* . One interesting question is the behavior 
of the system at low densities. 

This system attracted interest in both the theoretical physics and mathematics litera- 
ture, starting from the mid-1990's. It appears that the invariant distributions (Jl.lj) have 
been discovered several times independently. Among the early ones was Evans Q] who 
derived the invariant distributions for the disordered exclusion model in both continuous 
and discrete time. Independently, Krug and Ferrari [7] studied the phase transition of 
the continuous-time model and interpreted the results in various physical contexts such as 
traffic flow and directed polymers. In general, on the physics side there is wide interest in 
particle systems as simple models of traffic flow and other "single file" systems. We refer the 
reader to Nagel [S] for a review of particle systems in traffic modeling. The state-of-the-art 
in traffic modeling with exclusion type systems is the Gray-Griffeath model which is 
an exclusion process whose jump rates depend on nearby sites. 

Returning to the disordered exclusion, on the mathematical side Benjamini et al. [5j first 
proved hydrodynamic limits for several asymmetric exclusion and zero-range processes with 
random rates. However, their assumptions specifically ruled out the phase transition. 

A complete hydrodynamic limit theorem for the model studied here was proved by 
Seppalainen and Krug For the case p* > the result was the following. If the 

initial distributions have a macroscopic profile below p* , then on the hydrodynamic scale 
the initial macroscopic profile is rigidly translated at speed c. In particular, if the system 
has initially a spatially homogeneous particle distribution with density p < p* (such as 
ergodic gaps with mean u > u*), a tagged particle satisfies 

t~ 1 o'i(t) — > c as t — > oo. 

Subsequently Andjel et al. proved a weak convergence result for the low density regime. 
Start the system so that the gaps are ergodic with mean u > u* . Then the gap process 
converges weakly to the maximal invariant distribution, in other words to the product 
distribution with marginals as in with a = c. 

The hydrodynamic limit and the weak limit suggest the following picture. Let us follow 
particle <jq that initially starts at the origin. The other particles are distributed so that the 
gaps are for example i.i.d. with mean u > u* , and then initially particle density is p < p* . As 
t grows, particle ao(t) experiences an increasing density around itself, and correspondingly 
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its advance is slowed down. The reason is that <jq is part of an ever-growing "platoon" 
of particles, headed by an especially slow particle. As this platoon catches up with slower 
platoons ahead of it, it grows and slows down even more. As t — > oo, the particle density 
around <Jo(i) approaches the critical density p* , and simultaneously his motion slows down 
to rate c. However, all this must happen at a scale below the hydrodynamic, because the 
hydrodynamic limit reveals only the trivial final behavior. 

The purpose of this paper is to quantify the slowdown experienced by <7o(i) when the 
system starts at low density. Technically speaking we are seeking the next order term 
in the hydrodynamic limit. We find that by time t, <Jo(t) has traveled a distance ct + 
w(t)t( l>+1 ^ 1 / ( I/+2 ) where v > is an exponent characterizing the tail of F(p) as p \ c, and 
w(t) is a random quantity, which becomes strictly positive and is tight as t — ► oo. We do 
not have a precise limiting distribution for w(t). Our bounds suggest that for large t the 
tail of w(t) behaves like expj — C(u — u*)~ l w 2+u } for some constant C . These results are 
for annealed distributions, in other words for probabilities where the random rates have 
been averaged out. 

Following the nonrigorous picture sketched above, proofs of the estimates proceed by 
bounding the rate of the slowest particle in a suitable range ahead of &o(t). The technical 
side of the proofs involves couplings of various kinds between several processes with different 
rates and/or initial distributions. 

We also address another question which is related, and partly uses the same tools for the 
proof. When the exclusion process starts with all sites in (— oo,0] occupied and all sites in 
[1, oo) vacant, on the hydrodynamic scale there is a limiting macroscopic profile supported 
on ( — oo,c]. In other words, the number Xt of particles that are in (ct, oo) at time t must 
satisfy Xt = o(t). We find bounds on the true size of Xt. This question is not restricted to 
the situation where u* < oo. It makes sense whenever F{c) = because then every particle 
is attempting to jump at a rate strictly higher than c. Then presumably X t is unbounded 
as t increases. 



disordered exclusion process 
2. The results 
The basic assumption is on the tail of F(p) as p \ c. 

There exist constants — 1 < v < oo and < k < oo such that 
(2 ' 1 ' lim FW - 



p\c (p - C)^ 1 

If the reader prefers a concrete example, let F have density f(p) = k(v + l)(p — cf on 
some interval (c, c + e). At f = — 1 the distribution i 7 has a jump of size k at c, so there 
is a positive density k of particles with minimal rate c. The behaviors we look at become 
simple. Values v < — 1 are of course not possible. Recall the definition (|1.2[) of the critical 
gap u*. An integration by parts checks that, under assumption (|2.1|) . > is equivalent 
to u* < oo. 

First we look at the slowdown phenomenon in low density. We specify that particle 

o"o starts at the origin [ctq = 0]. Initial locations (<7j : i ^ 0) of the other particles are 

determined by taking the initial gaps {r^} i.i.d. random variables with common mean u = 

Erji > u* and finite variance. Then set 

i-l -l 
o"j = ^2 Vj f° r * > 0, and cr, = r/j for i < 0. 

j=0 j=i 

Our results are bounds on the "annealed" distributions of the quantities of interest. This 
means that while the process is run with fixed rates p = {pi}, we look at the average of 
all the processes for different choices of p, but with the fixed initial distribution for (o"i). 
P will denote this probability measure which represents the random choice of rates, the 
random initial configuration (<7j), and the random exclusion evolution. 
Notationally it is convenient to use 

1 



a 



u + 2' 

so that in particular the power of the correction is 

1 — a = . 

v + 2 

Set also 

A(v) = 'r—rrr. 
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Theorem 1. Assume (|2.1|) with v > 0. Let the initial gaps {r}{\ be i.i.d. random variables 
with common mean u = Era > u* and finite variance. The following bounds are valid for 
any < z < oo . 

m ^ P p(^ > ,) < «p{-^)-'^ *«} 

and 

(2.3) lisLnfP (-|^ >2 )> exp {__^_ 2 -}. 

Next we consider the situation where initially all sites in (— oo,0] are occupied by parti- 
cles, and all sites in [1, oo) are vacant. This could be thought of as an outflow from a large 
jam. Now there is always a rightmost particle, so we label the particles with nonpositive 
integers in increasing order. We drop the generic a notation, and for this special situation 
denote the locations of the particles at time t by 

■••<£_ 2 (t)<£-i(t) <£<>(*)■ 

The initial locations are £i(0) = i for i < 0. Particle £j jumps at rate pi independently 
drawn from distribution F. 

This system has a hydrodynamic limit which can be expressed in terms of the empirical 
measure as follows: for a compactly supported continuous test function <p, 



lim V^(t _1 ^(t)) = / 4>{x)r{x)dx 
t^roo — ' In 
i<0 Jtl 



almost surely. The limiting density r(x) is supported on (— oo,c]. (The reader can find 
more information about the limit and r(x) in ^J.) For the homogeneous exclusion with 
constant rates 1 this is Rost's classical result 9 , with a piecewise linear profile 



r\{x] 



1, x<-l 
5(1-2;), -Kx<l 



0, x > 1. 

The random rates produce the following qualitative difference with the homogeneous 
case. In the homogeneous case the lead particle £o(t) is a Poisson process of rate 1, and 
so its location is t + O^ 1 / 2 ). In other words, its location coincides with the right edge of 
the hydrodynamic front. However, in the disordered system the lead particle is a Poisson 
process of rate po, which under assumption ()2.1|) is strictly greater than c. Thus £o(t) and 
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in fact a large number of particles are ahead of the hydro dynamic front whose right edge at 
time t is at ct. The second question we address is to bound the number of these particles. 
Let X t be the number of particles that are beyond point ct at time t, in other words 

X t = 0V sup{k > 1 : £_fc+i(t) > ct}. 

Theorem 2. Assume (|2.1j) with v > 0. Then for all b > 0, 

limsupP{X t > bt l - a } < ex.v{-A{v)- l Kb u+2 \ 

and 

liminf > bt l - a } > expf-A(u)(l + u*) v+1 Kb v+2 \ 

When i/ < we no longer have a finite critical gap size u*. Theorem ^ fails, not just 
because u > u* is no longer possible, but because in equilibrium o" (t) is a Poisson process 
and has fluctuations on the scale i 1 / 2 . 

The phenomenon described by TheoremEJis not restricted to v > 0. With —1 < v < it 
is still the case that many particles advance ahead of the hydrodynamic front, as no particle 
has the lower bound c as its actual rate. 

For v = our result is the same as for v > but with a logarithmic weakening in the 
lower bound. This seems an artifact of our proof, so it is not clear whether this is the true 
state of affairs. Note that at v = we have a = 1 — a = 1/2, matching with diffusive 
fluctuations. 

Theorem 3. Assume (|2.1[) with v = 0. Let e > be arbitrarily small and < a < oo 
arbitrarily large. If < b < oo is large enough, then for all large enough t, 

P{at 1/2 {\ogty l <X t < bt 1/2 } > 1 - e. 

Xt changes behavior for u < 0, and is of smaller order than 0(t l ~ a ). Unfortunately we 
do not have matching upper and lower bounds. As v \ — 1 (a /* 1) the ratio of the upper 
and lower bound exponents becomes one. 

Theorem 4. Assume (12.1(1 with —1 < v < 0. Let e > 0. // < b < oo is large enough, 
then for all large enough t, 

P{X t < 6t (1+iy)/2 } > 1 - e. 
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// < a < oo is small enough, then for all large enough t, 

P{X t >at^/^}>l-e. 

The upper bound is on the boundary of conflicting with Gaussian fluctuations of the 
Poisson clocks. For large b, with high probability the slowest particle among bt^)l 2 
particles has rate at most c + qt~ l l 2 for a small q > 0. Consequently the number of jump 
attempts experienced by this slow particle by time t is Poisson with mean ct + qt 1 / 2 . This 
can be brought below ct by a fluctuation of order t 1 / 2 in the clock. Thus there is some 
chance that this particle does not reach ct by time t. To improve the probability to 1 — e we 
choose b and q so that there is a large enough number of slow particles. The lower bound 
meets this "Gaussian border" only in the limit v \ — 1. 

3. Variational representations 

In this section we run through notions which have been elaborated elsewhere flT. . The 
purpose is to establish the conventions followed in this paper which in some cases deviate 
slightly from those used before. Let an arbitrary initial configuration a = {cxj} be given, 
random or deterministic. Fix the rates {p{\. The process a(t) = {o~i(t)} is constructed with 
the usual graphical representation, by attaching a rate pt homogeneous Poisson process 
iVj = (Ni(t) : t > 0) to each particle <7j. 

Construct an auxiliary family {C(t)} of exclusion processes by stipulating that at time 
t = their initial locations are 

Cj(0) = a i +j forj<0. 

Only particle indices j < are used for the auxiliary processes. The jumps of the particles 
Cj are defined by 

Q attempts to jump whenever Poisson clock Ni+j rings. 

This translation of the index of the clock has the effect that for any fixed k, particles 
{°~ki Cjfc-j : * — ^} m ake jump attempts at the same times, namely when clock N)~ rings. 

Process has initially all sites in (— oo, Uj] occupied and all sites in [<7j + 1, oo) vacant. 
From this observation one can see that the variational equation 

(3-1) *k(t) = .inf d-M 

i:i>K 
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is valid at t = 0. Then one proves it by induction on jumps for all times t. 

In Theorems |2HH we consider the system £(t) that starts exactly as C(t) but centered at 
the origin. Let 

Then the processes £ % (t) are copies of except that the rates {pi} have been shifted in 
space. Of course this does not affect the distribution of when the rates are averaged 
out. We will find it convenient to use the variational equality ()3.1|) also in the form 

(3.2) a k (t)= infte + SU*)}- 

i:i>k 

Exclusion processes can be represented by interface processes. Suppose an interface 
process is given in terms of a height function i h,i{t) from Z into Z. This means that at 
time t the interface is the graph of the function h(t), so that hi(t) is the vertical coordinate 
of the location of the interface over site i. We impose the condition hi < hi+i on admissible 
height functions. Dynamics are defined by stipulating that if N{(t) = Ni(t—) + 1, then 

hi(t) = hi(t—) + 1, provided hi(t—) < hi+i(t—) — 1. 

In other words, height hi jumps up at rate pi, provided it does not go above its right 
neighbor. Obviously, we can map between a(t) and h(t) by 

Gi(t) = hi(t) + i. 

Precisely speaking, if the processes a{t) and h(t) are coupled so that this equality is true 
at t = 0, then it remains true for all t > 0. 

The gap process rj(t) = {r]i(t)} is defined in terms of these processes by 

Vi(t) = &i(t) ~ 0i-i(t) - 1 = hi(t) - hi-i(t). 

The variational equation for the height process takes this form. Let Z l (t) be an interface 
process with these properties: initially 

Zj = for j < i, and Zj = oo for j > i. 

Dynamically, 

Zj takes its jump commands from Poisson process Nj, for all i and j. 
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Then 

(3.3) h k (t) = inf {hi + Zlit)}. 

i:i>k 

There is no translation in ()3.3j) because each column of the height processes h(t) and Z l (t) 
reads the same clock. Since <Jo(t) = ^o(*)> we can use the variational formula 

(3.4) a (t)= wf {hi + Zitt)} 

i:i>0 

in the proof of Theorem ^ where we follow the evolution of o~o(t). 

4. Proof of Theorem [T] 

We begin with the key lemma that points the way to controlling the behavior of the 
system by looking at the slowest rate in a suitable range of indices. For fixed positive qi 
and q 2 , and a positive real parameter N, let 

(4.1) J(N) = >0: Pi <c + q 2 N- a } 
and define the event 

(4.2) D(N) = \pi > c+ q 2 N~ a for < i < [qiN 1 ^ = {J(N) > giiV 1 ^} . 

Lemma 1. Assume (|2.1j) and recall that the rates {pi} are i.i.d. with common distribution 
F. For fixed qi,q 2 > 0, 

(4.3) lim P[D(N)) = expj-^!^ 1 } . 



Proof. Let 5 > 0. For p sufficiently close to c, 

(4.4) (k -5)(p- cY +1 < F(p) <(k + S)(p - c) v+l 
Due to the independence of the rates pi, we have 

(4.5) P(D(N)) = {l-F(c + q 2 N- a )) [9lNl ' a] . 
This yields the upper and lower bounds 

1 - (« ± 5)q v 2 +l N- a{ - v+ ^ N 



for P(D(N)). Let N oo and then 5 to obtain the limit (fOjl . □ 
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4.1. Proof of the upper bound in Theorem^ The upper bound (|2.2j) follows from 
this proposition. 

Proposition 1. Suppose the initial gaps {r/.;} are an i.i.d. sequence with common mean 
u > u* and finite variance. Let qi,q 2 > 0. Then for any 5 > 0, 

lim sup P[a (t) >ct + ( qi (u - u*) + q 2 )t 1 ' a + St 1 ' ] < lim P(D(t)). 



Before proving Proposition^ let us observe how it implies the upper bound Q2.2|) . To- 
gether with Lemma n the proposition gives 



lim sup P [cj (t) > ct + zt 1 ' } < exp{-Kgig2 +1 } 

t— >oo 



for any q\, q% such that z = q\{u — u*)+q2 + S. Minimize the right-hand side of the inequality 
subject to this constraint on qi,q 2 - Then let S — > 0. 
The remainder of this section proves Proposition^ 

Lemma 2. Consider an arbitrary process o~(t). Let K > 1. Then 

1. 



lim P 



a (t) = min {hj + Zl(t)\ 
o<j<Kt J u 



Proof. From the definition of the process 2T- ^(-)> initially at time zero z\ (0) = for 
i < [ift]. Variable Z^-q is the first to jump, after which Z^j_ 1 may jump, then zjj^*-j_ 2 , 
and so on. Consequently the time T when variable Z^^ takes its first jump up is a sum of 
independent exponential waiting times with rates p\jcth P\Kt\-l-> P[Kt]-2i ■ ■ ■ > Po- Let e > 0. 
Since each rate pi is bounded above by one, T < (K — e)t with probability that vanishes 
exponentially fast as t — > oo. If we take < e < if — 1, we conclude that 

p{zf\t) >0}^0 

exponentially fast. 

To prove the lemma, it suffices to show that Zq ' (i) = implies 

<Jo{t) = min [/&,- + Zl(t)\. 
o<j<Kt J u 
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Let i > [Kt]. Then, since Zq > always and the height hi is nondecr easing in i, 

hi + 4(t) >hi> h [Kt] = h m + Zf l \t). 

This shows that indices i > [Kt] cannot contribute to the infimum in the variational 
formula. □ 

Lemma 3. Consider two processes a and a whose initial gaps are i.i.d. with common 
mean Era = Efji = u and finite variances. Couple the initial configurations so that they are 
independent, but give the processes the same rates {pi} and the same Poisson clocks. Then 
for 5 > 

lim P[a (t) > o o (t) + 5t l ~ a ] = 0. 

t^>oo 



Proof. Let K > 1 and define the event 



A(K, t) = <T (t) = mm + Z> (t)} and a Q (t) = mmlhj + Z J (t)} 

( 0<j<.Kt 0<j<Kt 

By Lemma! P(A(K,t)) -> 1 as t -> oo. On A(K,t), 

an(t) = min {h» + Zl[tY\ = min {Ha — hi + Ha + Zl(t)\ 

0<j<Kt J u 0<j<Kt J J J u 

< er (t) + max {hj r- hj}. 

0<3<Kt 



By Kolmogorov's inequality, 
P 



max -f/t,,- — hA > <5i 



l-a 



< 



Kt Varfai] + .Frt Var[r/i] 



^2(2(1-0) 

As 2(1 — a) = 2(f + l)/(f + 2) > 1, this last expression vanishes as i — > oo. Consequently 

P[(r (t) > er (t) + St 1 ' } < P(A(K, t) c ) + P 
gives the conclusion by letting t — » oo. □ 



max {/i,- — /io} > 5t 

0<j<Kt 



l-a 



Now define a particular mean u initial system as follows. Fix a number u < u* , and let 
a be the equilibrium velocity corresponding to average gap u, defined by 

dF(p). 



( c ,i] P - a 
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For each realization p of the rates, let {rji} have the nonstationary geometric product 
equilibrium distribution 



Then Er/i = u, and the {rji} are i.i.d. when the random rates are averaged out. We chose 
u strictly less than u* because then 



This finite variance is necessary so we can apply the previous Lemma |3 We cannot use 
equilibrium gaps at mean u* because they have infinite variance if < v < 1. 

Let {ji} be an i.i.d. sequence of nonnegative integer valued random variables, indepen- 
dent of {fji}, and with common mean Eji = u — u. Assume the {"fi} have finite variance. 
Define 



Then {fji} are i.i.d. with common mean u. Let a(t) denote the process with cto(O) = and 
initial gaps {fji}. 

Lemma 4. For any 5 > 0, 



Proof. Thinking of the zero-range process of the gap evolution, couple the processes fj(t) = 
{fji(t)} and fj(t) = {fji(t)} via the basic coupling, so that fji(t) < fji(t) for all i and t. This 
entails having cjj and (Xj read the same Poisson clocks for each i. 





fji =Vi + 7i- 



limsupP[a (t) > ct + (qi{u - u) + q 2 )t 1 ~ a + St 1 ""} < limsupP(£>(*)). 




Let 



J(t) = inf{» > : pi < c + q 2 t~ a }. 



J(t) depends only on the rates. Since ajr t \(t) — oj(i)(0) is stochastically dominated by a 
mean ct + q2t 1 ~ a Poisson random variable, the event 



Bi(t) = {a J{t) {t) < cr J(t) (0) + ct + q 2 t 1 ~ a + 5^/4} 



satisfies P(^i(t) c ) -» 0. Let 



Bz(t) = {a J{t) (0) < J(t)(u + 1) + J{t)8l{Aq x )}. 



14 ILIE GRIGORESCU, MIN KANG, AND TIMO SEPPALAINEN 

By the weak law of large numbers, P(B2(t) c ) — > because J(t) — > oo almost surely. By the 
connection between particles cij(t) and gaps and by the coupling with fji(t), 

J(t)-i 

ao(t) = aj(t)(*)- E ^(iW(i)<<7/(t)(*)- fji{t)-J(t). 

i=0 i=0 

By stationarity, = {^(i)} has the same distribution for all t > 0, under any fixed p. 
Now combine the inequalities. On the event 

A{t) = {or (t) > ct + (gi (u - u) + ^i 1 "" + Jt 1 ""}, 

we have 

j(t)-i 

*j(t)(*)> S ^W + ^)+ct + (gi(n-u) + g 2 )t 1 - Q + ( 5t 1 - a . 

i=0 

Consequently on n -Bi(i) we have 

*j (t) (o)> ^ ^) + J(t) + ?i(«-«)t 1 - a + ^t 1 - Q . 

j=0 

Next, on A(t) n n B 2 (*) we have 

Y, Vi(t) < J(t)u - qi (u- u)t l ~ a - -St 1 -* + 
i=o qi 

And finally, on the event D(t) c , J(t) < git 1- ", and so as our last inequality, on A(t) n 
Bi(t) n B 2 (t) n D(t) c we have 

J(t)-i 

m(t)<J(t)u--8t 1 -<*. 



2 



To summarize, 



P[a (t) >ct+ ( qi (u -u) + q 2 )t 1 - a + St 1 " } < P(D(t)) + P(B 1 (tY 
+P(B 2 (ty)+p(D(t) c n{ ]T J(t)« -^-72}) ■ 



i=0 



The conclusion follows because on D(t)°, J(t) < q\t a while still J(t) — ► 00, so the last 
probability vanishes as t — > 00. □ 
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Now we prove Proposition ^ Fix u < u* so that 

qt(u* - u) < 5 /4. 

Define the processes a(t) and a(t) as was done for Lemma HJ Couple all three processes 
(cr(t),a(t),a(t)) so that the initial gaps of o~(t) are independent of the initial gaps of the 
other two, and all read the same Poisson clocks. By the choice of u, 

P[° {t) >Ct+ ( qi ( U - U*) + q 2 )t 1 ~ a + (ft 1 ""] 

< P[*o(t) >ct + ( qi (u -u) + q 2 )t 1 - a + 5t l ~ a /A] 
+ P[<j Q {t)>d {t) + 8t l - a /A}. 

Let t — ► oo and apply the lemmas. 

4.2. Proof of the lower bound in Theorem^ The lower bound will follow from proving 
this proposition. 

Proposition 2. Suppose the initial gaps are i.i.d. with common mean u > u* and finite 
variance. Given positive q\,q 2 let 

r = min{(/2, (u — u*)qi}. 

Then for any 5 > 0, 

(4.6) liminf P{a Q (t) > ct + rt x ~ a - St 1 ' 01 } > expj-^i^+H . 



The lower bound (j2.3|) will follow from this proposition the same way the upper bound 
()2.2j) followed from Proposition ^ Namely, for a given z maximize the right hand side of 
()4.6j) subject to r — 5 = z, and then let 5^0. 

To prove Proposition [21 we start with the variational equation and split it into two 
separate ranges. 

a (t) = inf{^(0) + Zj»(i)} = mm{S 1 (t),S 2 (t)} 

where 

Si(t)= inf {hj(0) + Z J Jt)} and S 2 (t) = inf {hj(0) + Z 3 (t)} . 
0<i<git 1_Q j>git 1 ~ a 
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We shall show that 

(4.7) liminf P{SUt) > ct + q 2 t l ~ a - 8t l ~ a } > lim P(D(t)) 

t—KX> t^QO 

and 

(4.8) lim P{S 2 (t) >ct + qi (u- u*)t l ' a - St 1 '*} = 1. 
Together with Lemma ^ these imply Q4.HJI . 

Proof of lower bound, part 1. In this section we prove (j4.7|) for S\{t). 

Proposition 3. Let qi,q 2 ,5 > 0. There exists an event B{t) such that P(B(t) c ) — > and 

{S^t) >ct + q 2 t 1 - a - dt 1 - } D D(t) n B(t). 

Lower bound H4.7|) follows from this proposition. The rest of this subsection proves the 
proposition. Pick a further constant (73 such that 

< q 3 < Q2 < Q3 + 5/4. 
We shall couple o~(t) with a faster process a(t) whose jump rates pi are given by 

Pi= Pi V (c + q 2 N- a ). 

Process a(t) will be in equilibrium so that each particle &i(t) jumps as a Poisson process 
with rate 

a = c + q 3 N- a . 

To achieve this, the gap process fj(t) = {fji(t)} has to have the appropriate geometric 
product equilibrium distribution. Given p, {fji} are independent with geometric marginals 

i-B, = *l = (i- •)(»)', , = 0.1,2,... 

Note that this is sensible because a < pi for each i by the assumption q 3 < q 2 . The processes 
fj{t) and a(t) depend on N, but we suppress this dependence from the notation. 
The mean gap for the &(t) process is 

u = E[r)i] = [ ^-rdF(p). 

J(c,l] P~ a 

Lemma 5. The mean gap u converges to u* as N — > 00. 
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Proof. The integral comes in two parts: 



u 



I J^ dF{p)+ j ^dF{p) 

J{c,c+q 2 N- a ] V — a J(c+q 2 N- a ,l] V ~ a 

c + q 3 N- a KT-^.f c + q 3 N- a n , pn 
F(c + q 2 N )+/ — l( c+ g 2iV — ,i]W dF (P)- 



(q 2 - q 3 )N 01 -' y (C)1] p-c- q 3 N- 

The first term on the last line vanishes as N — > oo by hypothesis (|2.1|) . To the second 
term we apply dominated convergence. The integrand converges to c/{p — c) for each fixed 
p £ (c, 1] , and satisfies the bound 

c + q 3 N- a q 2 c + q 3 

: 1 (c+q 2 N-",l}{P) S 



p-c - q 3 N 01 { 92 ' J q 2 - q 3 p-c 

if N > 1. The last upper bound is integrable under dF(p), again by assumption (|2.1|) , □ 

For higher moments of fji we develop a bound. 

Lemma 6. For k > 1 and iV > 4, 

E[rft] < CN^- 1 "^ log N. 

C is a constant that depends on k and all the other constants in the problem, but not on 
N. 

Proof. For a fixed p, properties of a geometric distribution give 

E p [r)h < Co + C^m" < Co + d(pi - d)- k 
for constants Cq,C\ that depend on fc. It remains to show 

(p - a)- fc dF(p) < CiV Q(fe ~ 1 - v)+ log iV. 



'(c,l] 

This integral is decomposed as 

(4.9) (q 2 - q 3 )- k N ka [ dF(p) + [ (p - a)- k dF(p). 

J{c,c+q 2 N- a \ J(c+q 2 N- a ,l] 

Apply assumption Q2.1JI to the first integral. In the second integral observe that 

/ \ k 

(p ~ d)- k < q 4 (p - c)' k for g 4 = ( 

.Q2-q 3 , 

Subsume the constants qi into constants Cj. Thus the next upper bound is of the form 

(4.10) CoN^- 1 "^ + Ci / (p- c)~ k dF[p). 

J{c+q 2 N- a ,l] 
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Pick C2 and 5 > so that F(jp) < C2(p — c) u+l for c < p < c + 5. In the second term, the 
integral over (c + S, 1] is bounded by a constant. Over (c + q%N~ a , c + 5] integrate by parts. 

/ ^%<F(c + 6)S- k - f F(p)d{(p-c)- k } 

J(c+q 2 N- a ,c+S] \P - c ) J(c+q 2 N~ a ,c+5} 

<c 3 + c 2 k [ c+s {p- c y- k dp. 

Consider different cases for the last integral. If v > k — 1, it is bounded by a constant. If 
v = k — 1, it is bounded by C3 + C4 log N. Finally in the case v < k — 1 it is bounded by 
CsN^- 1 -"^. In all cases the bound given in the statement of the lemma works. □ 

Couple {rji} and so that they are mutually independent. 

Lemma 7. For any q > 0, 5 > 0, 

lim Pi inf [hi - L] < -SN 1 -"} = 0. 

iVote t/iat the height function h changes with N in the statement above. 

Proof. Take N large enough so that u — u > 0, which can be achieved by Lemma[21and the 
assumption u > u* . Then the probability in the statement of the Lemma is bounded above 
by 

Pi inf [h j -h j -j(u-uj\<-5N l - a \ 

< 5- 2 N- 2 ^- a ^ ■ 9 jV 1-a (Var[7fo] + Var[%]) 

where we used Kolmogorov's inequality. By the previous lemma Var[r/ ] < CN a ^-^ + logiV 
while Var[?7o] is a constant. As a(l — v) + < 1 — a for all v > 0, the probability vanishes as 
N -> 00. □ 

Now we turn to 5i(t). Consider first a fixed Set N = t, and as above construct the 
equilibrium process &(■) with rates p\. Also, let Z J denote the corner processes run with 
the pi rates. On the event D(t), we have 

Si(t) = min {h j + Zi(t)}= min {h j - L + L + Ztit)} 
0<j<gitl-« O^Xqii 1 "" J J 7 u 

> min {/i 7 - — /i 7 | + <7o(i) 

O^Xqit 1 "" 
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because 

hj + M{t) > & (t) 

for each j > 0. Consequently 

{5i(t) > ct + q 2 t 1 ~ a - 5t l - a ] 
3 Dft) n | min {hi - L} > -±<ft 1_a , <7 (*) > ct + g 2 i 1_a - iff* 

= D(t)nB(t), 

where the last identity means that the event B(t) is defined by the previous expression in 
braces. For the complement 

P(B(tr)<p{^_ a {^--/i i }<-^--} 

+ P {a (t) <ct + q 2 t 1 ~ a - ISt 1 -"} . 

The probabilities above vanish as t — > oo, the first by Lemma [7| above. For the second 
probability, note that <5"o(t) is Poisson distributed with mean 

at = ct + qst 1 ' >ct + q 2 t 1 ' a - \5t l ~ a . 

Since 1 — a > 1/2, the deviation j5t l ~ a has zero probability in the t — * oo limit. This 
completes the proof of Proposition |31 

Proof of lower bound, part 2. In this section we complete the proof of Proposition |^1 by 
proving (|4.8|) . 

Proposition 4. Given e, 8 > 0, 

P{S 2 (t) >ct+(u 

for all large enough t. 
Proof. Let {rj*} be the independent mean u* equilibrium gaps, so given p, 

p-k^ (!_£)(£)', k >_ . 

Let <r*(t) be the equilibrium process where particle c>o(i) is a rate c Poisson process. Couple 
the processes a(t) and <r*(i) so that they read the same Poisson clocks but their initial states 
are independent. 



- u* - 5)q 1 t 1 - a } > 1 - e 



20 ILIE GRIGORESCU, MIN KANG, AND TIMO SEPPALAlNEN 

By the Strong Law of Large Numbers, 

M _1 (hM — h* M ) > u — u* almost surely. 

Note that here we do not need finite variance, which the h* height function would not 
possess if < v < 1. Shrink 5 if necessary so that < 5 < u — u*. Pick Mq = Mq(5, e) such 
that 

P{hj - h* > j(u - u* - 5/2) for all j > M } > 1 - e/2. 
Since 1 — a > 1/2, there exists a io such that 

P{°o(t) >ct- q 1 t 1 - a 5/2} > 1 - e/2 

for all t > i . Now with probability at least 1 — e, for t > to such that qit 1-01 > Mq, 

S 2 (t) = inf {h j + Z i {t)} = inf {hj - h* + K + Z>(t)} 

> inf {h j - h*} + a*Jt) >q 1 t 1 - a (u-u*-5/2)+ct-q 1 t 1 - a 5/2 
j> qi ti-°< J 



ct + qi t L - a (u-u* -5). 



□ 



5. Proof of Theorem E] 

We begin with the upper bound. Let b> 0, O<0<1, 0<e< Ob and q 2 = Ob — e. Let 
p be the minimal rate among p_r w i-ai, . . . , P-je^i-a], and / an index such that pi = p. Let 
Y(t) be a Poisson variable with mean ct + q 2 t 1 ~ a ''■ If P < c + q 2 t~~ a , Y(t) dominates the 
number of jump attempts particle experiences during time interval [0,t]. By the particle 
ordering, ^r^i-ai (t) > ct implies > ct, and thereby £1 must have at least ct + Obt 1 ^ 01 
jump attempts. We get the bound 

P{X t > bt l - a ] < P{^ [u i- a] (t) > ct} 

<P{p>c + q 2 t~ a } + P{Y(t) >ct + 0bt l ~ a }. 

Since 1 — a > 1/2, the last probability vanishes as t — > oo. By Lemma ^ we get 

limsupPLXt > bt l ~ a ] < exp(-«(l - a\u^+^ 



i2 

t— >oo 



Let e \ so that q 2 f Ob, and then choose = (v + l)/{v + 2). 
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The lower bound of Theorem Incomes from the lower bound of Theorem^ Pick a density 
u > u* , and let the initial gaps {rji} be bounded i.i.d. random variables with mean u. By 
the variational formula (|ci.2|) . 

a (t) = M{a j + e(t)} 

where is a version of the process with translated rates. Let b > 0, and then pick 
9 > b(u + 1). Let j = [bt l - a \. Then 

gLfi) > a (t) - Oj >ct + (a (t) - ct) - Uj. 

The annealed distribution of the process £ J (t) is the same as that of £(t). Consequently 

p{x t > bt 1 -*} > p{^ [btl - a] (t) > ct} > p| ao( t l; ct > e, a lbt i- a] < et 1 -"}. 

By the law of large numbers t +Q! <7r(, t i-c«i — > b(u + 1), and so by Theorem^ 

liminf P{X t > bt l - a \ > liminfp{ a ° ( ?~ C * > e] > expl —8» +2 \. 

Maximize the last lower bound over 9 and u subject to u > u* and 9 > b(u + 1). 

6. Proof of Theorem El 

The argument for the upper bound is similar to the previous one. Now v = and 
1 — a = 1/2. Let e > be small. By the central limit theorem we can fix a large 
1 < M < oo so that, if Y(t) is a Poisson random variable with mean ct + t 1 ~ a , then 

P[Y(t) > ct + Mt l - a ] < e/4 

for all large enough t. Given e and M, choose 0<q2<l<M<q<bso that 

exp(-K^ +1 ) > 1 - e/16 and exp(-K&<^ +1 ) < e/16. 

Let 

J(t) = inf{i >0:p-i<c+ q 2 t~ a }. 
By LemmaElwe have to < 00 so that 

P^ 1 -" < J(t) < to 1 " } > 1 - e/4 

for all t > to- Suppose this event happens. Then if £_[^i-a](t) > ct, also £_j( t )(t) > 
ct, and particle C-J(t) nas had to cover distance ct + J(t) > ct + qt l ~ a . The increment 
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£-j(t)(t) — £-j(t)(0) is stochastically bounded by the variable Y(t) defined above. So for 
large enough t, 

P{X t > bt 1 -*} < P{£- [bt i- a] (t) > ct} < P{£_ J(t) (i) > ct} + e/4 

< P{Y{t) >ct + qt l - a ] + e/4 < e/2. 

We prove the lower bound by comparison with a faster system in equilibrium. Let 
< a < oo be fixed. Given e > 0, pick 1 < w < oo large enough so that 

(6.1) P[Y(N) > EY(N) - wN 1/2 } > 1 -e/4 

for large enough N, for a Poisson variable Y(N) with mean cN + 2WN 1 / 2 . Later we have 
to increase w further. 

Let q2 = 4io, and define faster rates by pi = pi V (c + q2N _1 / 2 ). Consider iV large enough 
to have pj < 1. Let a(t) be a process run with rates pi and in equilibrium, so that &o(t) is 
a Poisson process with rate 

a = c + 2wN~ l/2 . 

The gap process fj(t) then has a product distribution with independent geometric marginals 

P- = = (!-£)(£)', , = 0,1,2,... 
The annealed mean gap is 

u = E[fji] = [ dF(p) 
J(c,i] P-a 

and the annealed variance is bounded as in 

Varft] < E[fjf] = 2 / (^—) dF(p) + u. 

J(c,l] \P - a J 

Lemma 8. There is a constant C that depends only on the distribution F such that for 
large enough N , 

u<ClogN and Varfa-] < CN 1/2 . 

Proof. First for the mean. Integrate by parts, and use assumption (j2.1j) to pick < 5 < 1 
such that F(p) < (re + l)(p — c) for c < p < c + 5. Then note that 

p — a > — (p — c) = h(p — c) for p > c + q2N~ 1 / 2 . 

Q2 
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Carrying out these steps yields 

aF(c + q 2 N- 1 / 2 ) f a 
u = zz 1 /a \- 



I dF{p) 
Jfc+ao.N- 1 / 2 ,!] P~ a 



c + q 2 N- 1 / 2 -a J(c+q 2 N- 1 / 2 ,i] P ~ a 
aF(c + q 2 N~ 1 / 2 ) 



C + Q2N- 1 / 2 



a 



I- a c + q 2 N~ l l 2 - a J( c +q 2 N-^,i} \P ~ a 



a 



a f 1 F(p) , 
+ a / ^dp 



r+() dp f F(p) 



a f c+ ° dp f 

<- + 4(/e + l)o/ — — +4a/ 

1 - a J c +g 2 iv-i/2 p-c J c+S [p - c) 2 

< — h 4(k + l)a flog (5 - log g 2 iV _1/2N ) + 4a<5~ 1 

1 — a V / 

< + 2(k + 1) log N + 45' 1 

< ClogN. 

In the second last step we took N large enough so that 

a = c + q 2 N~ 1/2 < — ^ < 1. 
If TV > 3, in the last step we can take 

1 + C . _ . , x . 



C 



— - + 2(k + 1) + 4<T 



which depends only on the distribution F. 

Following the same pattern for E[fjf] shows that, after integration by parts, the main 
part is the integral 

" C+S dp , o 2 f 1 F(p) 



2 f 1 Hp) , 2( , n r +5 ^ ^ 8 2 r 

a / - dp < 8a (k + 1) / 7 + 8a / 

Jc+qiN-^ [P ~ a) d J c+q2N -i/2 (p - cy J c+S 



(p - c) d 



The desired bound follows as above. □ 
Let denote a £-type process run with rates j5j. Let 



j(N) = [aN 1 / 2 (log N)' 1 }. 
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From the variational coupling (|H.'2j) we have 

= ct + 2wN~ 1/2 t + (a (t) - (ct + 2wN~ 1/2 t)) - a j{N) . 

The processes and a(t) depend on A" but we suppress this from the notation. Set time 
t = N. Note that when the random rates are averaged out, processes |*W(t) &nd£(t) have 
the same distribution. We get this bound. 

P{L j{N )(N) > cN} 

> p{a (N) - (cN + 2wN 1 / 2 ) > -wN 1 ' 2 , a j{N) < wN 1 ' 2 } 

> P{a (N) - (cN + 2wN 1/2 ) > -wN 1/2 } - P{a j{N) > wN 1/2 }. 

The next to last probability is at least 1 — e/4 for large N by . It remains to show that 
the last probability vanishes as N — » oo. From the annealed perspective <?j(N) ls a sum °f 
i.i.d.'s, so its mean and variance are bounded as follows: 

Ea j{N) = j(N)(u + 1) < CaN 1/2 and Vax[a j{N) ] = j(N) Var[r/ ] < CaN(log N)' 1 . 

At this point we need to increase our original choice of w to guarantee that w > 2Ca where 
a is given in the beginning of the proof and C is the constant that appears in Lemma |H1 
Then Chebychev's inequality gives 

P{* m > W NV 2 } < P{a j{N) > Ea j{N) + CaN^ 2 } < 

which vanishes as — > oo. We can conclude that for large N, 

P{L j{N )(N) > cN} >l-e/3. 

Finally we make contact with ^(A r ). Given qi chosen above, pick q\ > small enough so 
that exp(— Kq\q2) > 1 — e/7. Let D(N) be the event 

D(N) = { Pi = pi for —[qiN 1 / 2 ] <i<0}. 

By LemmalU P(D(N)) > 1 - e/6 for large enough N. On the event D(N), ^(t) = |j(t) 
for —[qiN 1 / 2 ] < i < and all t > 0, so in particular for i = j(N) if A^ is large enough. 
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Consequently 

P{X N > aN 1 / 2 (log N)- 1 } > P{£_j (N) (N) > cN} 

> P({t- J(N) (N) > cN} n D(N)) = P[{i- J(N) (N) > cN} n D(N)) 

> P{€-j(N)(N) > cN} - P(D(N) C ) > 1 - e/3 - e/6 = 1 - e/2. 
This completes the proof of Theorem |3J 

7. Proof of Theorem |U 

7.1. Proof of the upper bound of Theorem |31 The upper bound is proved by com- 
parison with independent particles. Let 

I — a 1 + 2/ 

/3 = ^r = — 

For b > and #2 > 0, define 

o 

if t = Y, HPi<c + q 2 t- 1/2 }. 
i=-[btf]+i 

Lemma 9. Let {Yj(t)} be independent copies of a Poisson random variable with mean 
ct + q2t 1 / 2 , independent of the rates {pi} and thereby independent of Kt. Then given e > 0, 
if Q2 is small enough while bq^ +l is large enough, 

P{Yj{t)>ct forl<j <K t ) <e 

for all large enough t. 

Proof. Fix a small < 5 < 1/2. Fix a positive integer m large enough so that {\+o~) m < e/2. 
Pick Eq > small enough so that 

P(X > so) < (1 + S)/2 

for a standard normal x- Let qi < Eq^Jc. 

By assumption (|2.1jl . for large t Kt is stochastically dominated by a binomial random 
variable with [bt 13 ] trials and success probability (k + l)q% +1 t~ ( 1+I/ )/ 2 . Such a variable 
converges weakly to a Poisson with mean b(n + 1)^2 +1 as * ~* 00 • Thus we may fix 6 large 
enough so that 

P(-Kt < "») < e/2 

for large enough i. 
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By the choice of q 2 and the definition of Y(t) = Yi(t), 



1 K ' - \ VVaxy(t) 



Then by the central limit theorem, for large enough t 



P{Y(t) > ct} < P{ X > so) + 5/2 < 1/2 + S. 



Finally, as the Yj(t) are i.i.d. and independent of K t , 



K t 



P{Yj(t) > ct for 1 < j < K t } = E Yl p {Yj{t) > ct} 

Lj=i 

< E\{\ + 5) Kt ] < P{K t < m) + (I + <5) m < e. 



Fix b and ^2 so that the lemma is satisfied. Let 



h = {-[bt 13 ] <i<0: Pi < c + q 2 t- 1/2 } 



□ 



Once the rates pi have been chosen according to distribution F and It determined, give 
each index i G It an independent Poisson process TVj(-) of rate c + q 2 t~~ 1/2 . Thin ATj(-) 
appropriately to get the correct rate pi. These thinned processes are the Poisson clocks for 
indices i G It- Meanwhile, give the other indices their independent Poisson clocks. This way 
we can claim that for each i G It, the number of jump attempts experienced by particle 
£i during (0, t] is bounded above by the mean ct + q2t 1 ^ 2 Poisson variable iVj(i) that is 
independent of the rates pi. 

Suppose £__[£,£/3] + i(i) > ct. By the particle ordering, £i(t) > ct for all i G It, which implies 
that Ni(t) > ct for all i £ It- By the lemma above this event has probability less than e for 
large t. To summarize, we have shown that for an arbitrary e > 0, b can be chosen so that 

P{X t > bt p } < e 



for large enough t. 
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7.2. Proof of the lower bound of Theorem |1J For an exclusion process with constant 
rates r, for any a > and < 7 < 1, 

(7.1) Hm ^^7^ = -2VH in probability. 

This statement is a consequence of a limit proved by Glynn and Whitt |Hj and the explicit 
computation of the value 2 on the right hand side first done in JU]. See Lemma 4.1 in |1U| 
for the derivation of (|7.1jl from [H]. (But note that the process £ in ^Oj is not the same as 
£ in the present paper.) 

Let j3 = (3 + ^) • Let < a < oo and (7 = 2-^/a + 2. Use assumption (|2.1j) exactly as in 
the proof of Lemma^fc show that, given e > 0, if a is small enough, then for large enough 
t 

P{Pi > c + qt'P for -[at /3(1+iy) ] < i < 0} > 1 - e/2. 

On this event ^_r a ^(i+f)](t) is bounded below by £_[at/3( 1 + i 0](*) where is a process whose 
clocks ring at constant rate c + qt~@. For ()7.1() gives the following bound: for large t 
with probability at least 1 — e/2, 

L [at M+>»](f) > <* + 1^ ~ 2 Mct + qt 1 -?)^ 1 ^ 2 ~ t^ 1+ ^ 2 > Ct. 

The last lower bound by ct followed from 1 — f3 = (1 + /3(1 + ^))/2 and the choice of g. 

We have shown that, given s > and a small enough a > 0, then for large enough £, the 
inequality C-[ a t/ 8 ( 1 + i/ )](*) > c ^ holds with probability at least 1 — e. This inequality implies 
X 4 >at^ 1+ ^. 
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